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Abstract 

We derive a (Wilsonian) effective field theory of QCD at finite density by 
integrating out the states in the Dirac sea when the chemical potential /j, 3> 
Aqcd- The quark-gluon coupling is effectively (l+l)-dimensional and the 
theory contains four-quark operators which become relevant as we approach to 
the Fermi sea. By calculating the one-loop vacuum polarization tensor in the 
effective theory, we find the electric gluons have a screening mass, M ~ g s fi, 
while the static magnetic gluons are unscreened. We then investigate the gap 
equations for color anti-triplet Cooper pairs by including both gluon-exchange 
interactions and the marginal four-quark interactions in the effective theory. 
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Quantum chromodynamics (QCD) has become a undisputed theory for the strong inter- 
action of quarks and gluons. However, its direct experimental evidences are still limited to 
the high energy region where quarks and gluons are weakly coupled. At low energy where 
they couple strongly, we have only indirect evidences due to the calculational incapability 
that QCD describes the strong interaction. 

Recently, QCD at high density has been studied intensively, since it not only offers a 
nonperturbative test but also exhibits rich novel phases which might be accessible in the core 
of neutron stars or in a quark-gluon plasma in heavy ion collisions [l|-|6| • It has been argued 
that at a sufficiently high density the color symmetry is spontaneously broken, leading to 
color superconductivity. 

In this letter, we derive a (Wilsonian) effective theory of QCD at high density in loop 
expansion, after integrating out the irrelevant degrees of freedom at low energy. The ef- 
fective theory includes new quark-gluon couplings generated by the exchange of states in 
the Dirac sea, which induce marginal four-fermion operators at matching one-loop ampli- 
tudes, as in the effective theory of QED under external magnetic fields 0. The marginal 
four-quark operators become relevant as we scale down to the Fermi sea. In the effective 
theory, the gluon-quark couplings are effectively (l+l)-dimensional and the electric gluons, 
Aq(x), are screened in the static limit. At energy below the screening mass, the only rel- 
evant interactions for quarks are the four-quark operators with opposite momenta and the 
interactions with magnetic gluons. We then solve the gap equations for Cooper pairs in the 
color anti-triplet channel by long-range color magnetic interactions, screened color electric 
interactions, and the marginal four-quark interactions. 

A system of degenerate quarks with a fixed baryon number is described by the QCD 
Lagrangian density with a chemical potential /i, 

£ QCD = ^P^- \F; u F a ^ + ^7o^, (1) 

where the co variant derivative = + iggA^T 11 and we neglect the mass of quarks for 
simplicity. (Later, we will consider the effect of quark mass.) 

For a sufficiently high density such that \i ^> Aqcd, one can derive systematically a 
(Wilsonian) effective theory of QCD, where Aqcd/a* provides a useful expansion parameter. 
At an energy scale just below the chemical potential /i (3> Aq CD ), quarks are almost free 
and their energy spectrum is given by an eigenvalue equation, 

(a-p- fj,)ip± = E ± ip ± , (2) 

where a = 707 and ip± denote the energy eigenfunctions with eigenvalues E± = —fi ± \p\, 
respectively. At low energy E < fi, the states ip + near the Fermi surface, \p\ ~ /x, are easily 
excited but ■?/>-, which correspond to the states in the Dirac completely decoupled 

due to the presence of the energy gap /j provided by the Fermi sea. Therefore the right 
degrees of freedom below \i consist of gluons and ip+ only. 

Since at very low energy (E -C /i) quarks are almost on-shell and move at the Fermi mo- 
mentum pf = fiVF, we may decompose the momentum of quarks into the Fermi momentum 
and a residual momentum as 

P^ = fiv^ + In, (3) 
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where the residual momentum < (i and t> M = (0,v F ). Now, we decompose the quark 
fields as 

= £ [e^-^ + {v F , x) + e^-^VF, x)] , (4) 

Vp 

where a ■ v F il)±(v F ,x) = ±ip±(v F ,x). The quark Lagrangian in Eq. (|Xj) then becomes 
4) (i p + Ai7°) ^ = E y j +(^F,x)i^D fJi ilj + (v F ,x) + i>-(v F ,x)-f° [2/j, + iD\\J ip-(v F ,x) 

+ E [$-(v F ,x)i P±Mv F ,x) +h.c. (5) 



where 7 f = (7°, ^ • 7), 7? = 7 M - 7f > ^11 = F"D„ with = (1, -tfr), and £ ± = t££> m . 

At low energy, since the fast modes ^_ are decoupled, we integrate out all the fast modes 
ip- to derive the low energy effective Lagrangian by matching all the one-light particle 
irreducible amplitudes containing gluons and tp + in loop expansion. The effects of fast 
modes will appear in the quantum corrections to the couplings of low energy interactions. 
At tree-level, the matching is equivalent to eliminating tp_ in terms of equations of motion; 

ry° ry° 00 / iD\\\ n 

i[)-(v F ,x) = -- — — — P ± ij + (v F ,x) = -— ^ ~o P±^+(v F ,x). (6) 
2{jt + iD\\ 2/i 2(i J 

Therefore, the tree-level Lagrangian for ip + becomes, using (1 — a ■ v F )/2^{l + a ■ v F )/2 = 
7°V*(1 + a ■ v F )/2 with V = (1, v F ), 



C = E 



(7) 



where the ellipsis denotes terms with higher derivatives. We see that in the high density 
limit the quark-gluon couplings are effectively (l + l)-dimensional: In the leading order, to 
quarks of Fermi velocity v F only the V ■ A component of gluons couple and = (0, A) 
with v F • A — does not couple to the color charge of ip + but to its color magnetic moment, 
which is suppressed by 1///. 

By the tree-level matching shown in Fig. 1, we obtain a vertex of two gluons and two 
quarks, generated by the exchange of ip-, 

££ = "f|Etf+AL7oAiV4- + -", (8) 

" Vp 

where the ellipses denote terms containing more powers of gluons and derivatives. (From 
now on ip will denote the slow modes ip+.) 

Now, consider the one-loop matching of a four-quark amplitude. (See Fig. 2.) In QCD, 
the amplitude is ultraviolet (UV) finite but infrared (IR) divergent, while it is both UV and 
IR divergent in the effective theory. Since the IR divergence is same in both theories, we 
need a UV counter term in the effective theory to match the amplitude, which is a four-quark 
operator. The one-loop four-quark amplitude in the effective theory is 
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where / is the external momentum transfer, A is a renormalization point, and c.t. denotes 
the counter terms. u,s,t,v are color indices and S^ s . tv = (5 uv 5 ts + 5 ut 5 sv ) /y/2, 5^ s . tv = 
(SuvSts — S u t5sv) I V%- We find that the new quark-gluon coupling generates effective four- 
quark operators at one-loop; 
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At the matching scale A = /x, the couplings are of order of a 2 . As in the BCS superconduc- 
tivity P|, this four-quark operator is marginal under the scale transformation toward the 
Fermi surface, lu — > s/ii with s < 1, only when the incoming quarks have opposite Fermi 
momenta, p\ + P2 = + 0(1), where we write the quark momenta pf = /iff + /f as before. 
Therefore, if this four-quark operator leads to a condensate, it must be in the s-channel, 
i'4 ; s{Pi)' l Pu{P2)) 7^ or (ip s (vp,x)ip u (—vp,x)) ^ 0. Note also that we do not have marginal 
operators containing ^(vp, x)tp s (vp, x) because it vanishes identically for the slow modes 
ip+. Therefore, one can easily understand in the effective theory that the chiral symmetry 
will be restored as the chemical potential becomes large, if any translationally invariant 
order parameter for chiral symmetry breaking should contain the fast modes ■?/>_ ||. 

Next, we match the gluon two-point amplitude at one-loop. Since the quarks in the Dirac 
sea are integrated out, we need to match only the quark- loop correction to the amplitude. 
The quark-loop contribution to the gluon two-point function in QCD at finite density and 
at zero temperature is calculated by Manuel [fLOR . It consists of two parts; one due to the 
matter and the other due to the vacuum. In terms of the vacuum polarization tensor, 



n^,fuii(p) 



where n^ vac is the quark-loop contribution when there 



is no matter, \i = 0. For \p^\ <C /i, the matter part of the vacuum polarization becomes for 
N f light quarks, with M 2 = N f g 2 fi 2 /(2n 2 ), 
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which is transversal, p M n^ mat (p) = fig. 

On the other hand, in the effective theory the quark-loop correction to the vacuum 
polarization is given as: 



1 — a ■ Vp 
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To match the gluon two-point amplitudes in both theories, we therefore need to add a term 
in the one-loop effective Lagrangian, 

(13) 

Vp 

which also ensures the gauge invariance of the effective Lagrangian at one-loop. Since in 
the effective theory II°° ~ —iM 2 5 a b, as p — > 0, we see that the electric gluons, A%, have a 
screening mass, but the static magnetic gluons are not screened at one-loop due to the added 
term Eq. ([T3|), which holds at all orders in perturbation as in the finite temperature |TT| , p"2l| . 

Finally, we match the quark two-point amplitudes to get a one-loop low energy (Wilso- 
nian) effective Lagrangian density 

1 2 A/f 2 - 1 

C eS = --(1 + en) - —AfA a ^ + (1 + h)^D^ - — (1 + (l± • S) 2 V 
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(finite - rf S ;tJ (#f> %)iPs(v F , x)tpl(-v F , x)ip u (-v F , x) (14) 



+ (hSt;tv ~ h e S ns;tv) V4(#F, x)^ a (v F , x)<lpl(-V F , x)j 5 <lp u (-V F , x) 



+ • • 



where the summation over v F is suppressed and the coefficients ai,&i,Ci are dimensionless 
and of order a s (fi). The ellipsis denotes the irrelevant four-quark operators and terms with 
more external fields and derivatives. 

As we scale further down, the effective four-quark operators will evolve together with 
other operators, which can be seen by further integrating out the high frequency modes, 
Sfi < {Ipl < ft. The scale dependence of the four-quark operators has three pieces. One 
is from the one-loop matching condition for the four-quark amplitudes and the other two 
are from the loop corrections to the four-quark operators, shown in Fig. 3. Putting all 
contribution together, we find the one-loop renormalization group equations for the four- 
quark operators to be 

where i = (6,3), g 6 = -g 6 , = g% and 7; = (a/2/9) (13/4, 1/2) and 84 = (-1,2). Since 
in the high density limit the quark-gluon coupling is (l+l)-dimensional, the quarks do not 
contribute to the running of the strong coupling. The one-loop (3 function for the strong 
coupling constant at high density is /3(a s ) = — ll/(27r)a 2 . Integrating the RG equation (|15|), 
we find for fi ^> A 



2ir 

&(A) ~ — a a (A) 



7i 144 1 



(16) 



At a scale much less than the chemical potential, A <C ft, (73(A) ~ 0.04a s (A) and g&(A) ~ 
0.29a s (A). Similarly for h h ^ = \/2/2(2, 1) and S { = (-1,2) and we get hg(A) ~ 0.4a s (A) 
and h 6 (A) ~ 0.81a fl (A). 

At a scale below the screening mass, we further integrate out the electric gluons, which 
will generate four-quark interactions. For quarks moving with opposite Fermi momenta, the 
electric-gluon exchange four-quark interaction is given as 
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A 9 3 -^^Y,^° T Mvf,x)^ TM-Vf,x). (17) 
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Using Tt u T£ s = l/28 ts S uv — l/65 tu 5vs, we find that the four-quark couplings are shifted 
as g- 3 {M) -> g- 3 (M) + 2 v / 2£ 2 (M)/3 ~ 0.95(? 2 (M) and # 6 (M) -> <? 6 (M) + - 
0.49g 2 (M), while hi's are unchanged. 

As we approach further to the Fermi surface, closer than the screening mass M, the four- 
quark operators in the color anti-triplet channel become stronger, because the (3 function 
for the attractive four-quark operators is negative, /3(g 3 ) = — g 2 /(47r 2 ). If the four-quark 
interaction is dominant at low energy, it leads to vacuum instability in the infrared region 
by forming a color anti-triplet condensate or Cooper pair, e tsu /^"'(up, x)ip^(— vp, x)\ ^ 0, 
where flavor indices a, (3 and Dirac indices i,j are restored. The size of the condensate 
driven by the four-quark interaction is determined by the RG invariant scale, 

( 4vr 2 \ 

A-Mexp 57T7r • (18) 

But, since the long-range color-magnetic interactions become also strong at low energy, we 
need to consider both interactions to determine the Cooper-pair gap. 

To describe the Cooper-pair gap equation, we introduce a charge conjugate field, 

{i>c) a = C a(} Mx), (19) 

where a and (3 are Dirac indices and the matrix C satisfies C _1 7 M C = — 7J. Then, we may 
write the inverse propagator for *ff(v F ,x) = (■ip(vF,x),ipc(—VF,%)) T as 

s-^-o^Uw t. ( #)- < 20 > 

The gap equation for A in the color-antitriplet channel is then given in the hard dense loop 
(HDL) approximation as (see Fig. 4) 

A(P||) " { ~ %9s) J j2^ D ^ iP ~ l)V / 2 -A 2 (Z„) V + V J (2 7 r)M;-A 2 (/ || )' (21) 

where D^ u is the gluon propagator. 

The gluon propagator is given in the HDL approximation as,[| following the notations 



used by Schafer and Wilczek |L5 



' d ^-Sg + Sf-^ < 22 > 



where £ is the gauge parameter and the projectors are defined by 



In the Schwinger-Dyson equation the loop momentum should take the whole range up to the 
ultraviolet cutoff, which is the chemical potential /i in the case of high density effective theory. 
Hence the gluon propagator includes both magnetic and electric gluons. 
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In the weak coupling limit, \k§\ <C \k\ and thus 
F(k ,k)~M 2 , G(k ,k) 



~ —M 2 -S- 

4 \k\ 



(25) 



Since the gap has to be fully antisymmetric in color indices, we get 



T t a u A uv (T a ) 
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After Wick-rotating into Euclidean space, the gap equation becomes 



(26) 



A(P||): 
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Note that the main contribution Q to the integration comes from the loop momenta in the 
region gjj ~ A 2 and \q±\ ~ M 2 / 3 A 1,/3 . Therefore, we find that the leading contribution is by 
the first term due to the Landau-damped magnetic gluons. For this momentum range, we 
can take \p — q\ ~ \q±\ and 



V ■ P T ■ V 
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J\>ikj 
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(28) 



We also note that the term due to the four-Fermi operator is negligible, since g% ~ g 4 at the 
matching scale \x. 

Neglecting (p — q) 2 in the denominator, the gap equation becomes at the leading order 
in the weak coupling expansion and l//i expansion 



m) 3 J (2tt) 4 



+ 



qt + W 2 \Po-qo\/\q±\ q± + M 2 
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The q± integration can now be performed easily to get 



An Mm) 

(2ir) 2 q 2 + A- 
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(30) 



2 The hard dense loop approximation is therefore consistent since A <C M. I thank Steve Hsu for 
explaining this point. 
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We see that in this approximation A(pii) ~ A(po)- Then, we can integrate over vp • q to get 
A(ni ) = -^r*b-^Inf * 1 (31) 



36tt 2 J-p ^ q 2 + A 2 V |p - go 
where A = 4///7T • (/i/M) 5 e 3/2? . If we take A ~ A(0) for a rough estimate of the gap 
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or A^Aexpl- — J. (32) 



As was done by Son [|T^], one can convert the Schwinger- Dyson gap equation ([Tl]) into a 
differential equation to take into account the energy dependence of the gap. Approximating 
the logarithm in the gap equation as 



we get 



hi bo Qo\ j i n |g | otherwise, 



pA-W + A'W + I^^O, (34) 



where p = po- The solution to the differential equation (33) is found to be at p = A(p) 



7T 



A~Aexp(^-- + l + 0(z/ 2 )J, (35) 

where z/ = \J%a s / (97r). The gap is therefore given as at the leading order in the weak 
coupling expansion?] 

A=c i exp B6)' (36) 

where c = 2 7 tt 4 NJ 5 ^ 2 e 3 ^ 2 ^ +1 . This agrees with the RG analysis done by Son |TJ|] (see also ||14|| ) 
and also with the Schwinger- Dyson approach in full QCD lfL6| , 15|, 17|| . The l/g s behavior of 
the exponent of the gap at high density is due to the double logarithmic divergence in the 
gap equation (Wn), similarly to the case of chiral symmetry breaking under external magnetic 
fields |I8| , |T9|| . In addition to the usual logarithmic divergence in the quark propagator as in 
the BCS superconductivity, there is another logarithmic divergence due to the long-range 
gluon exchange interaction, which occurs when the gluon loop momentum is colinear to the 
incoming quark momentum (q± — > 0). 

The color anti-triplet condensate has an interesting flavor and Lorentz structure as dis- 
cussed in [fj|^,[12| . So far we have neglected quark masses, since they are irrelevant at 



3 The gauge-parameter dependent term is sub leading in the gap equation (|27]) . Since the gap has 
to be gauge-independent, the gauge parameter dependence in the prefactor will disappear if one 
includes the higher order corrections. 
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scales higher than the quark masses, namely at A > m s , the mass of strange quark in the 
three flavor case. But, if g§ or g s do not become strong enough to form a condensate at 
the scale A = m s , the strange quark becomes irrelevant in the low energy effective theory 
and decouples from the condensate. Furthermore, if the condensate does not form at low 
energy scale where the instanton effects are important, one may expect the instanton contri- 
bution to the condensate will be quite significant fHH. But, as analyzed in || in detail, the 
instanton generated operators are not dominant when the chemical potential is quite large. 

In conclusion, we have derived a (Wilsonian) low energy effective theory of QCD at high 
density in loop expansion by integrating out the antiparticles. The effective theory contains 
marginal four-quark operators which become relevant as we approach to the Fermi sea. We 
also calculate the screening mass for gluons and find that electric gluons have screening 
mass, M ~ g s fi, while the magnetic gluons are not screened in the static limit. Finally, 
in the effective theory, we solve the gap equation for the Cooper pairs in the color anti- 
triplet channel in the hard dense loop approximation. We find that the gap is given as 
A ~ fig~ 5 exp[— 3n 2 /(\/2g s )] up to a numerical constant in the prefactor, verifying recent 
results. 

ACKNOWLEDGMENTS 

The author is grateful to Sekhar Chivukula, Andy Cohen, Nick Evans, Steve Hsu, Roman 
Jackiw, Rob Pisarski, Dirk Rischke, and Raman Sundrum for useful discussions. This work 
was supported in part by the KOSEF through SRC program of SNU-CTP and also by the 
academic research fund of Ministry of Education, Republic of Korea, Project No. BSRI-98- 
2413. 



9 



REFERENCES 



[1] D. Bailin and A. Love, Phys. Rep. 107 (1984) 325. 

[2] R. Rapp, T. Schafer, E.V. Shuryak, and M. Velkovsky, Phys. Rev. Lett. 81 (1998) 53. 
[3] M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B 422 (1998) 247, [Eep"- 



ph/9711395| ; J. Berges and K. Rajagopal, Nucl. Phys. B 538 (1999) 215, [Kep- 
ph/9804233| . 



[4] M. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B 537 (1999) 443, [Eepj 
ph/9804403| . 



[5] N. Evans, S. Hsu, and M. Schwetz, Nucl. Phys. B 551 (1999) 275, |hep-ph/980844l 



Phys. Lett. B 449 (1999) 281, |hep-ph/98105T4l ; T. Schafer and F. Wilczek, Phys. Lett. 

B 450 (1999) 325, frep-ph/ 98 10509] 
[6] K. Langfeld and M. Rho, Quark condensation, Induced symmetry breaking and Color 

superconductivity at high density, |hep-ph/9811227 . 
[7] D. K. Hong, Phys. Rev. D 57 (1998) 3759, |hep-ph/970745^ . 

[8] J. Polchinski, in Proceedings of the 1992 TASI, eds. J. Harvey and J. Polchinski (World 
Scientific, Singapore 1993); G. Benfatto and G. Gallavotti, J. Stat. Phys. 59 (1990) 541; 
Phys. Rev. C42 (1990) 9967; R. Shankar, Physica A 177, (1991) 530; Rev. Mod. Phys. 
66 (1993) 129. 

[9] An interesting possibility for chiral symmetry breaking at high density by color-flavor 
locking for the three light flavors is suggested in |4j]. 
[10] C. Manuel, Phys. Rev. D 53 (1996) 5866; G. Alexanian and V.P Nair, Phys. Lett. B 
390 (1997) 370. 

[11] See, for example, J. Kapusta, Field Theory at Finite Temperature and Density (Cam- 
bridge University Press, Cambridge, England, 1994). 
[12] R. Pisarski and D. Rischke, Phys. Rev. Lett. 83 (1999) 37, |nucl-th/98llTU4 . 
[13] D.T. Son, Phys. Rev. D 59 (1999) 94019. 
[14] S. D. H. Hsu, and M. Schwetz, |hep-ph/99083Tq . 
[15] T. Schafer and F. Wilczek, |hep-ph/99065T2| . 

[16] D.K. Hong, V.A. Miransky, LA. Shovkovy, and L.C.R. Wijewardhana, |hep-ph/ 9906478 . 
[17] R. Pisarski and D. Rischke, nucl-th/990741; |nucl-th/9910056| . 

[18] V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, Phys. Rev. D 52 (1995) 4747; Nucl. 
Phys. B 462 (1995) 249; C. N. Leung, Y. J. Ng, and A. W. Ackley, Phys. Rev. D 54 
(1996) 4181; D.-S. Lee, C.N. Leung, and Y.J. Ng, Phys. Rev. D 55 (1997) 6504. 

[19] D. K. Hong, Y. Kim, and S.-J. Sin, Phys. Rev. D 54 (1996) 7879. 

[20] V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, Phys. Rev. Lett. 83 (1999) 1291, 
ihep-ph/9811079| . 



10 



FIGURES 




FIG. 1. The tree-level matching condition. Wiggly lines denote gluons; solid lines, states near 
the Fermi surface; and double solid line, states in the Dirac sea. 




FIG. 2. The one-loop matching condition for a four-quark amplitude. 
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